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An empirical s c Pa tight-binding (TB) model is applied to the investigation of electronic states in 
semiconductor quantum dots. A basis set of three p-orbitals at the anions and one s-orbital at the 
cations is chosen. Matrix elements up to the second nearest neighbors and the spin-orbit coupling 
are included in our TB-model. The parametrization is chosen so that the effective masses, the spin- 
orbit-splitting and the gap energy of the bulk CdSe and ZnSe are reproduced. Within this reduced 
s cPa TB-basis the valence (p-) bands are excellently reproduced and the conduction (s-) band is 
well reproduced close to the F-point, i.e. near to the band gap. In terms of this model much larger 
systems can be described than within a (more realistic) sp 3 s*-basis. The quantum dot is modelled 
by using the (bulk) TB-parameters for the particular material at those sites occupied by atoms of 
this material. Within this TB-model we study pyramidal-shaped CdSe quantum dots embedded in 
a ZnSe matrix and free spherical CdSe quantum dots (nanocrystals) . Strain-effects are included by 
using an appropriate model strain field. Within the TB-model, the strain-effects can be artifically 
switched off to investigate the infuence of strain on the bound electronic states and, in particular, 
their spatial orientation. The theoretical results for spherical nanocrystals are compared with data 
from tunneling spectroscopy and optical experiments. Furthermore the influence of the spin-orbit 
coupling is investigated. 

PACS numbers: 73.22.Dj, 68.65.Hb, 71.15.Ap 
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I. INTRODUCTION 

Semiconductor quantum dots^ 2 * (QDs) are of particu- 
lar interest, both concerning basic research and possible 
applications. QDs arc considered to be zero dimensional 
objects, i.e. systems confined in all three directions of 
space with a typical size of the magnitude of several 
nanometers. Therefore, these systems are realizations 
of "artificial atoms" whose form and size can be manip- 
ulated. Concerning basic research these nanostructures 
(QDs) are interesting, as the methods of quantum theory 
can be applied to systems on new scales and with new 
symmetries in between that of atoms or molecules and 
of macroscopic crystals. On the other side light emission 
and absorption just from the localized states in such de- 
vices may be important for optoelectronic applications 3,4 , 
quantum cryptography 5 and quantum computing 6 . 

Semiconductor QDs can be realized by means of 
metallic gates providing external (electrostatic) confine- 
ment potentials?, by means of selforganized cluster- 
ing of certain atoms in the Stranski-Krastanow (SK) 
growth mod o 8 : 9 : 10 or chemically by stopping the crys- 
tallographic growth using suitable surfactant materi- 
alsiiiSi 1 ^* 1 ^. Here we deal only with the latter two types 
of QDs. The QDs created in the SK growth mode emerge 
self- assembled or self-organized in the epitaxial growth 
process because of the preferential deposition of material 
in regions of intrinsic strain or along certain crystallo- 
graphic directions. In epitaxial growth of a semiconduc- 
tor material A on top of a semiconductor material B only 
one or a few monolayers of A material may be deposited 
homogeneously as a quasi two dimensional (2d) A-layer 
on top of the B-surface forming the so called wetting layer 
(WL). Under certain conditions and for certain materi- 



als further deposited A-atoms will not form a further 
homogeneous layer but they will cluster and form islands 
of A-material because this may lower the elastic energy 
due to the lattice mismatch of the A- and B-material. If 
one then stops the growth process, one has free A-QDs 
on top of an A-WL on the B-material. If one continues 
the epitaxial growth process with B-material, one obtains 
embedded quantum dots (EQDs), i.e. QDs of A-material 
on top of an A-WL embedded within B-material. 
The chemically realized QDs emerge by means of colloidal 
chemical synthesisii* 1 ^. Thereby the crystal growth of 
semiconductor material in the surrounding of soap-like 
films called surfactants is stopped when the surface is 
covered by a monolayer of surfactant material. Thus 
one obtains tiny crystallites of the nanometer size in all 
three directions of space, why these QDs are also called 
"nanocrystals" (NCs). The size and the shape of the 
grown NCs can be controlled by external parameters like 
growth time, temperature, concentration and surfactant 
material 1 ^* 1 *. Certain physical properties like the band 
gap (and thus the color) depend crucially on the size of 
the NCs. Typical diameters for both, EQDs and NCs, 
are between 3 and 30 nm, i.e. they contain between 10 3 
up to 10 5 atoms. Therefore, EQDs and NCs can be con- 
sidered to be a new, artficial kind of condensed matter 
in between molecules and solids. For the in SK-modus 
grown EQDs lens-shaped dots 1 ^, dome shaped and pyra- 
midal dots&i&ii, and also truncated conesiS have been 
found and considered. 

Of course the fundamental task is the calculation of 
the electronic properties of EQDs and NCs. But here 
one encounters the difficulty that these systems are much 
larger than conventional molecules and that the funda- 
mental symmetry of solid state physics, namely trans- 
lational invariance, is not fulfilled. Therefore, neither 
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the standard methods of theoretical chemistry nor the 
ones of solid state theory can immediately be applied 
to systems with up to 10 5 atoms. Conventional ab- 
initio methods of solid state theory based on density 
functional theory (DFT) and local density approxima- 
tion (LDA) would require supercell calculations. But the 
size of a supercell must be larger than the EQD or NC, 
and such large supercells are still beyond the possibil- 
ity of present day computational equipment. Therefore, 
only systems with up to a few hundred atoms can be 
investigated in the framework of the standard ab-initio 
DFT methodo 19 i 20 i 21 . Simple model studies based on the 
effective mass approximation 1 ^^ or a multi-band k ■ p- 
mode! 23-24 25 describe the QD by a confinement potential 
caused by the band offsets, for instance; they give qual- 
itative insights into the formation of bound (hole and 
electron) states, but they are too crude for quantitative, 
material specific results or predictions. More suitable 
for a microscopic description are empirical pseudopoten- 
tial methodsSiiSLSSiS as well as empirical tight-binding 

models 30,31,32,33,34,35,36,37,38,39,40,41^ The empirical pseu _ 

dopotential methods allow for a detailed variation of the 
wave functions on the atomic scale. This is certainly 
the most accurate description from a microscopic, atom- 
istic viewpoint, but it requires a large set of basis states. 
Within a TB-modcl some kind of coarse graining is made 
and one studies spatial variations only on inter-atomic 
scales and no longer within one unit cell. The advantage 
is that usually a small basis set is sufficient, which allows 
for the possibility to study larger systems. Furthermore 
the TB-model provides a simple physical picture in terms 
of the atomic orbitals and on-site and inter-site matrix 
elements between these orbitals. A cutoff after a few 
neighbor shells is usually justified for orbitals localized 
at the atomic sites. 

Semiempirical TB-models have been used already to 
describe " nearly" spherical InAs and CdSe NCs for which 
the dangling bonds at the surfaces are saturated by 
hydrogen 31 i 32 i 33 i 34 or organic liganda 37 ! 38 ! 39 . Also un- 
capped 4 ^ and capped 3 ^ pyramidal InAs QDs were ives- 
tigated by use of an empirical TB-model. In the latter 
work an sp 3 s*-basis was used leading to a ION x ION 
Hamiltonian matrix, where N is the number of atoms, 
with 33 independent parameters. In the present paper 
we apply a similar TB-model to II- VI nanostructures, 
namely CdSe EQDs embedded within ZnSe and spher- 
ical CdSe NCs. We show that a smaller TB-basis is 
sufficient, namely an s c p 3 -basis, i.e. 4 states per unit 
cell and spin direction. This requires only 8 indepen- 
dent parameters and, in principle, allows for the inves- 
tigation of larger nanostructures than were accessible in 
Ref. |35[ Strictly speaking, the s c p 3 -basis-set leads to a 
smaller matrix-dimension and also to a smaller number of 
nonzero matrix elements compared to sp 3 s* TB-modcl. 
So the s c Pa TB-model is numerically less demanding re- 
garding both memory requirements and computational 
time. For the bulk system the valence p-bands are ex- 
cellently reproduced and the conduction s-band is well 



CdSe ZnSe 



E g [eV] 1.74S 2.820144 

A so [eV] 0.4li2 0.4342 

m e 0.12*2 0.14?44 

71 3.3344 2.4544 

72 1.11^ 0.6144 

73 1.45^ 1.1144 
C12 [GPa] 46.344 50.6^ 
Cn [GPa] 66.^ 85.944 



TABLE I: Properties of the CdSe and ZnSe bandstructures. 
The lattice constants are given by 6.077 A and 5.668 A, 
respectively. E s denotes the band gap, A so the spin-orbit 
coupling and m e the effective electron mass. The Kohn- 
Luttinger-Parameters are 71,72 and 73. The G%j are the ele- 
ments of the elastic stiffness tensor. 



reproduced close to the r-point. Therefore, we expect 
that also for the QDs all the hole states and at least the 
lowest lying electron states (close to the gap) are well re- 
produced. We investigate, in particular, the influence of 
strain effects on the electronic structure. To examine the 
accuracy of our model we compare the results to other 
microscopic and macroscopic models. Furthermore TB- 
results obtained for CdSe-NCs are compared to experi- 
mental results, and very good agreement, for instance for 
the dependence of the energy gap on the NC-diameter, is 
obtained. This demonstrates that our TB-model with a 
reduced basis set is reliable and sufficient for the repro- 
duction of the most essential electronic properties of the 
nanostructures. 

This work is organized as follows. In Sec. Irflour TB- 
model is presented. The formalism how to obtain the 
TB-parameters and how to apply them to the description 
of EQDs and NCs is described. In Sec. IHII the inclusion 
of strain effects in our model is introduced. Results for 
the pyramidal CdSe EQDs are presented. For the spher- 
ical CdSe NCs the results and the comparison with the 
experimental data are presented in Sec. IIVI Section [V] 
contains a summary and a conclusion. 



II. THEORY 

A. TB-Model for bulk materials 

In this work we use a TB-Model with 8 basis states 
per unit cell. Such a model has been succesfully used for 
the investigation of optical properties in ZnSe-quantum 
wells 4 ^. For the description of the bulk semiconductor 
compounds CdSe and ZnSe we choose an s c p\ basis set. 
That implies that the set of basis states \is, a, tr, R) is 
given by four orbitals a = s,p x ,p y ,p z with spin er = ±i. 
One s-orbital at the cation (y = c) and three p-orbitals 
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FIG 
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1: Tight- binding band structures for CdSe (a) and ZnSe 



at the anion (y = a) site in each unit cell R are chosen. 
The TB matrix elements are given by 

E a<a >(R' - R)„,»> = {v',a',o'R!\H hvSk \v,a,a,R) . (1) 

The coupling of the basis orbitals is limited to nearest 
and next nearest neighbors. Following Ref. |4?1 the spin- 
orbit component of the bulk Hamiltonian H hulk couples 
only p orbitals at the same atom. With the two center 
approximation of Slater and Kostes^ we are left with 
only 8 independent matrix elements E a ^ a i(R' — R) v , v i . 
In k space, with the basis states \k,i/,a,a), the electronic 
properties of the pure bulk material are modelled by an 
8x8 matrix H bulk (fc) (for each k from the first Bril- 
louin zone). This matrix depends on the different TB- 
parameters E aa i{R! — R)„ ]V '. By analytical diagonal- 
ization for special fc directions, the electronic dispersion 
E n (k) is obtained as a function of the TB-parameters; 
here n is the band index. Equations for the different 
TB-parameters E aa i(R' — R) v . v i can now be deduced as 
a function of the Kohn-Luttinger-parameters (71,72,73), 
the energy gap -E gap , the effective electron mass m e and 
the spin-orbit-splitting A so . The zero level of the en- 
ergy scale is fixed to the valence-band maximum. The 
used material parameters for CdSe and ZnSe are given 
in Table |U The resulting numerical values for the differ- 
ent TB-parameters (obtained by optimizing them so that 
the resulting TB band-structure reproduces the param- 
eters given in Table are summarized in the Table |n] 
(with and without taking into account a site-diagonal 
parameter for the spin-orbit coupling). Within this ap- 
proach, the characteristic properties of the band struc- 
ture in the region of the T point are well reproduced, 
as can be seen from Fig. ^ which shows the TB-bands 
of bulk CdSe and ZnSe (using the TB-parameters with 
spin-orbit coupling). When comparing with band struc- 
ture results from the literature^, one sees that the three 
valence (p-) bands are excellently reproduced whereas the 
s-like conduction band is well reproduced only close to 
the T-point. This is understandable, because higher (un- 
occupied) conduction bands are neglected, and can be 
improved by taking into account more basis states per 
unit cell. But for a reproduction of the electronic prop- 
erties in the region near the T-point, which is important 



Material 


Parameter 


TB 


TB-NO SO 


ZnSe 


E xx (000) aa 


-1.7277 


-2.0413 




E ss (000) cc 


7.0462 


12.1223 




E sx (lll)ac 


1.1581 


0.2990 




e** (no) aa 


0.1044 


0.2185 




E xx (011) aa 


0.1874 


0.0732 




^(110) oo 


0.3143 


0.4285 




Ess (110) cc 


-0.3522 


-0.7752 




A 


0.1433 





CdSe 


E xx (000) aa 


-1.2738 


-1.7805 




Ess(000) cc 


3.6697 


10.8053 




Esx (m) ac 


1.1396 


0.4260 




E xx (110) ao 


0.0552 


0.2161 




E*» (0H) aa 


0.1738 


0.0129 




E*y(no) aa 


0.1512 


0.3120 




Ess (H0) cc 


-0.1608 


-0.7554 




A 


0.1367 






TABLE II: TB-parameters (in eV) with (TB) and without 
(TB-NO SO) spin-orbit coupling for ZnSe and CdSe, using 
the notation of Ref. 48. 



for a proper description of the optical properties of the 
semiconductor materials, the s c p„-TB-model is certainly 
sufficient and satisfactory. 



B. TB-Model for Embedded Quantum Dots and 
Nanocrystals 

Having determined suitable TB-parameters for the 
bulk materials (here CdSe and ZnSe) a EQD or NC can 
be modelled simply by using the TB-parameters of the 
bulk materials for those sites (or unit cells) occupied by 
atoms (or molecules) of this material. Concerning the 
on-site matrix elements this condition is unambiguous. 
Concerning the intersite matrix elements one also uses 
the bulk matrix elements, if the two sites are occupied 
by the same kind of material, but one has to use suitable 
averages of the bulk intersite matrix elements for matrix 
elements over interfaces between different material, i.e. 
if the two sites (or unit cells) are occupied by different 
atoms (or molecules) . Concerning the surfaces or bound- 
aries of the nanostructure there are different possibilities. 
One can use fixed boundary conditions, i.e. effectively 
use zero for the hopping matrix elements from a surface 
atom to its fictitious nearest neighbors, or (for the em- 
bedded QDs) one can use periodic boundary conditions 
to avoid any surface effects, which artificially arise from 
the finite cell size used for the EQD-modclling. For the 
NCs the best thing to do is a realistic, atomistic mod- 
elling of the organic ligands covering the NC-surface, as 
described in Refsi22*£2iiii. Within the restricted basis set 
thus selected the ansatz for an electronic eigenstate of 
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the EQD or NC is, of course, a linear combination of the 
atomic orbitals \u, a, er, R): 



/ m, 5 1 v. a, a, i?) 

/ j i/,a,(7,itl 5 5)/ 



(2) 



Here R denotes the unit cell, a the orbital type, a the 
spin and v an anion or cation. Then the Schrodinger 
equation leads to the following finite matrix eigenvalue 
problem: 



(u ,a \a , R'\H\v,a,(7, R)u £—Eu , , , b, 

\ 3 3 3 II'''/ b>.o>,<j,R i/,a,cr,R' 



(3) 

where -E is the energy eigenvalue. The shortcut nota- 
tion (v' ,a' ,a' , R'\H\v,a,a, R) — H^, m g is used in the 
following for the matrix elements with / = v',a',a' and 
m = v, a, a. 

The matrix elements for CdSe and ZnSe without strain 
are denoted by H°~, For these matrix elements the 

J IR'.mR. 

TB-parameters E aa '(R' — i?)„,„< of the bulk materials, 
determined in Sec. Ill Al are used. For the off-diagonal 
matrix elements over interfaces and the diagonal matrix 
elements of the selen atoms at the interface between dot 
and barrier, which can not unambiguously be referred to 
belong to the ZnSe or CdSe, respectively, we choose the 
mean value of the parameters for the two materials. Fur- 
thermore, a parameter for the valence-band offset AEy 
has to be included in the model. This means that for 
CdSe in a heterostructure, i.e. surrounded by a barrier 
ZnSe material, all diagonal matrix elements are shifted 
just by AEy compared to the bulk CdSe diagonal ma- 
trix elements. In the literature different values for AEy 
can be found, they vary in the range of 10 %-30 % of the 
band gap difference between CdSe and XnSr' We 
have performed calculations with valence-band offsets of 
AE V = 0.108 eV, AE V = 0.22 eV and AE V = 0.324 
eV, which corresponds to 10 %, 20 % and 30 % of the dif- 
ference of the band gaps. We find that these different 
choices for AEy shift the EQD energy gap E^ by less 
than 2%. This shows, that the results are not much af- 
fected by the specific choice of the valence-band offset 
AEy- Therefore, in the following, an intermediate value 
of AE V = 0.22 eV is chosen. 

Furthermore, in a heterostructure of two materials 
with different lattice constants, strain effects have to be 
included for a realistic description of the electronic states, 
because the distance between two CdSe unit cells and 
the bond angles are not the same as the corresponding 
equilibrium values in bulk CdSe. This means that the 
intersite TB matrix elements H^, m g in the EQD differ 
from the - matrix elements in the bulk material. 

lR\mR 

In general, a relation 



tt ttO 
n lR' ,mR n iR'm 



Rf( d R>-R> d R>-R) 



(4) 



has to be expected, where d|, - and d^,_j^ are the 
bond vectors between the atomic positions of the un- 



strained and strained material, respectively. The func- 
tion f(d°,d) describes, in general, the influence of the 
bond length and the bond angle on the intersite (hop- 
ping) matrix elements. For lack of a microscopic the- 
ory for the functional form we use as a simplified 



model assumption /(^g 



R-> d R'-R) ~ { d %-R'-R'-Rj 

With this d~ 2 ansatz, the interatomic matrix elements 
^ir' m,w witn " ^ are S iveri by 



f)/ d R>-R 



0, 



TT — R° 

n lR',mR IRJ m R 



'd°, - v 

R'-R 
d R'-R, 



(5) 



This corresponds to Harrison's^ d~ 2 rule, the validity 
of which has been demonstrated for II-VI-materials by 
Sapra et alm^. More sophisticated ways to treat the scal- 
ing of the interatomic matrix elements, e.g. by calculat- 
ing the dependence of energy bands on volume effects and 
different exponents for different orbitals, can be found in 
the literatureSi^iiSi. Furthermore the results of Bertho 
et al*! for the calculations of hydrostatic and uniaxial 
deformation potentials in case of ZnSe show that the d~ 2 
rule should be a reasonable approximation. Our model 
assumption for the function f(S > ,d) means that we ne- 
glect the influence of bond angle distortion. Though en- 
ergy shifts due to bond angle distortions have been found 
for InAs EQDs^, here the negligence of bond angle dis- 
tortion can be justified when exclusively taking into ac- 
count the coupling between s- and p-orbitals at nearest 
neighbor sites. Piezoelectric fields, which are usually con- 
sidered to be less important for the zinc blende structures 
realized in CdSe and ZnSe24, are also not taken into ac- 
count in our model. 

The problem is now reduced to the diagonalization of a 
finite but very large matrix. To calculate the eigenvalues 
of this matrix, in particular the bound electronic states 
in the QD, the folded spectrum method^ is applied to 
the eigenvalue problem of Eq. |SJ 



III. RESULTS FOR A PYRAMIDAL CdSe 
EMBEDDED QUANTUM DOT 

A. Geometry and Strain 

To model a CdSe QD embedded into a ZnSe barrier 
material we choose a finite (zinc blende) lattice within a 
box with fixed boundary conditions. Within this box we 
consider a CdSe WL of thickness la (lattice constant of 
the conventional unit cell, i.e. about two anion and two 
cation layers), and on top of this wetting layer there is 
a pyramidal QD with base length b and height h = 6/2. 
For the matrix elements corresponding to sites within 
the WL or the QD we choose the TB-values appropriate 
for CdSe, for all other sites within the box the ones for 
ZnSe. Figure |21 shows a schematic picture of this geome- 
try we use to model the EQD. We investigate EQDs with 
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12.5 nm 



ZnSe^ 




CdSe wetting layer 



FIG. 2: Schematic visualization of the pyramidal CdSe QD 
buried in a ZnSe matrix. The wetting layer has a thickness 
of one lattice constant (la) of bulk ZnSe. The pyramidal QD 
has a base length b of ten times the ZnSe lattice constant 
(b = 10 a). 



a base lengths b of 6 a, 8 a and 10 a, where a = 5.668 A 
is the lattice constant of the bulk ZnSe material. Cells 
with the dimensions of 18 a x 18 a x 15 a (38 880 atoms), 
20 a x 20 a x 16 a (51 200 atoms) and 22 a x 22 a x 17 a (65 
824 atoms) are used for the calculations. Figure [21 shows 
the EQD with a base length b of 10 a. Fixed boundary 
conditions are applied to avoid a dot-dot coupling in con- 
trast to periodic boundary conditions—. The total size of 
the cells is chosen so that the boundary conditions affect 
the energy gap of the EQD by less than 2 %. 

To consider strain effects in our model the knowledge 
of the strain tensor e is necessary. The strain tensor e is 
related to the strain dependend relative atomic positions 

d A ,_ A =(l + e)d£,_ M . (6) 

To appoint the strain tensor outside the EQD, the WL 
is treated as a quantum film. In the absence of a shear 
strain (e^ j ~ <5.;.j) for a coherently grown film, the strain 
components are given byS 



e± = e zz = - — e|| . (8) 

Here oc is the lattice constant of the unstrained film ma- 
terial and as denotes the parallel lattice constant of the 
substrate. In Table [I] the cubic elastic constants Cy of 
the bulk materials are given. The resulting strain profile 
for a line scan in z-direction outside the dot is shown in 
Fig. (a). In Ref. [23| Stier et al. considered a similar 
strain profile for an InAs/GaAs EQD. The lattice mis- 
match of approximately 7% in the InAs/GaAs system is 
nearly the same as for the CdSe/ZnSe system. So our 



calculated strain profile shows the same behavior as the 
profile in Ref. |23j for a line scan in z-direction outside the 
EQD. 

To obtain the strain profile inside the EQD we use a 
model strain profile, which shows a simi lar beh avior as 
the strain profiles which are given in Refs. l23l6fl for a line 
scan in z-direction through the tip of the pyramid. This 
model strain profile is displayed in Fig. 0(b). The shear 
components, e xy , e xz and e yz , can be neglected, at least 
away from the boundaries of the dot2S. 



B. Bound single particle states 

We have calculated the first five states for electrons 
and holes for three different EQD sizes. These calcula- 
tions are done with and without including strain effects. 
For the evaluations without strain we have chosen the 
exponent in Eq. JSJ to be zero. The energy spectrum 
obtained from these calculations is shown in Fig. 0] (a) 
without strain and in Fig. 0] (b) including strain effects. 
The states are labeled by e\ and h\ for electron and hole 
ground states, ei and hi for the first excited states, and 
so on. All energies are measured relative to the valence- 
band maximum of ZnSe. Figure 0] also shows the size 
dependence of the electron and hole energy levels. The 
energies are compared to the ground state energies for 
electrons and holes in the 1 a thick CdSe WL (WL ei and 
WL/,j, respectively), which is calculated separately for a 
coherently strained quantum film (i.e. the WL without 
the QD). As expected from a naive particle in a box pic- 
ture, the binding of electrons and holes becomes stronger 
in the EQD when the dot size is increased. The quan- 
tum confinement causes the number of bound states to 




5 10 15 20 5 10 15 20 



z [nm] z [nm] 

FIG. 3: Strain distribution in and around the embedded pyra- 
midal CdSe QD with a base length of b = 10 a. The WL at 
the base of the QD is 1 a thick. The whole structure is buried 
in a ZnSe matrix. Line scans along the [001] direction through 
the WL outside the dot (a) and inside the dot through the tip 
of the pyramid (b) are displayed. The diagonal elements of 
the strain tensor e are shown as solid (e zz ) and dashed-dotted 
lines (e xx = e yy ). 
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Pyramid base length b [nm] 

FIG. 4: Electron and hole energies for embedded pyramidal 
CdSe/ZnSe QDs on a WL of thickness la (roughly 4 mono- 
layers) for different base length of the dot without strain (a) 
and with strain effects taken into account (b). The ground 
state energies for electrons (WL ei ) and holes (WLjJ for the 
WL (of thickness la) alone are also displayed. 



decrease when the dot size is reduced. For the EQDs 
with a base length b = 8 a and b = 10 a the calucated 
hole states are well above the WL energy (WL^ ) . This 
is valid for the strain- unaffected and strained EQD. For 
the system with b = 6 a we obtain at least four bound hole 
states in both models. The energy splitting between the 
different states is only slightly influenced by the strain. 
Furthermore we see from Fig.0|that the number of bound 
electron states is influenced by the strain. For the sys- 
tem with a base length of b = 10 a we get at least three 
bound-electron states when we take strain effects into 
account (Fig. 0] (b)). Without strain effects at least 5 
bound states are found. So the confinement potential for 
the electrons is effectively reduced by the strain. 

The bound electron states ei and are energetically 
not degenerate even without strain. This arises from the 
Ci v symmetry of the system. Already from the geometry 
of the EQD-system it is clear that there is no (001) mir- 
ror plane. Furthermore, if one considers a (OOl)-plane 
with sites occupied by Se anions, the nearest neighbor 
(cation) planes in ±z-direction are not equivalent, as in 
the zinc blende structure the nearest neighbors above the 
plane are found in [lll]-direction and below the plane 
in [lll]-direction. So also for crystallographic reasons a 
(OOl)-plane is not a mirror plane. Finally, if one considers 
the base plane of the EQD (or the WL) to be this anion 
(OOl)-plane, there are different cations, namely Cd above 
and Zn below this plane. Therefore, the QD-system has 
reduced C^-symmetry. In theories based on continuum 
models, e.g. effective mass approximationsi^S^, the dis- 
cussed effects cannot be accounted for. These inter- 
facial effects also affect the one particle wave functions 



in the system. In Fig. [3] the isosurfaces for the squared 
electron wavefunctions |$i(r)| 2 arc displayed with and 
without strain, respectively. The light and dark isosur- 
face levels are selected as 0.1 and 0.5 of the maximum 
probability density, respectively. For both calculations, 
the lowest electron state e\ is an s-like state according 
to its nodal structure. The next two states ei and are 
p-like states. These states are oriented along the [110] 
and the [110] direction, respectively. Due to the different 
atomic structure along these directions we find a p-state 
splitting A° 2 = E e3 — E &2 for the unstrained EQD of 

about 0.43 meV. In conventional k ■ p models 2 ^^ an un- 
strained, square-based pyramidal EQD is modelled with 
a C^ v symmetry. In our microscopic model the resulting 
degeneracy is lifted and a splitting occurs as a conse- 
quence of the reduction of C^ v symmetry to a C2 V zinc 
blende symmetry. 

The strain splits the states e2 and further. Due to 
the different atomic structure, the strain profile within 
each plane (perpendicular to the growth z-direction) 
along the [110] and [110] direction is different. This 
effect contributes also to the anisotropy. Due to the 
fact, that the base is larger than the top, there is a 
gradient in the strain tensor between the top and the 
bottom of the pyramid. In the EQD, the cation neigh- 
bors above each anion are found in [111] direction while 
the cation neighbors below are found in [lll]-direction. 
Therefore, the cations along the [110] direction are sys- 
tematically more stressed than the cations along the [110] 
direction. In case of strain we find a p-state splitting of 
A|'™ n — 7.1 meV. Compared to the states e-i and e 3 
of the unstrained EQD, the two lumps of the light iso- 
surfaces are well separated. The states and reveal 
nodal planes along the [110] and [110] direction, respec- 
tively. 

The state for the strained dot is resonant with a WL- 
state, so the wave function is leaking into the WL. Also 
the wave function of the state e§ is localized at the base 
of the pyramid but clearly shows already a finite proba- 
bility density inside the WL. The states e± and e$ of the 
unstrained EQD are still mainly localized inside the dot. 
The classification of the state en by its nodal structure is 
difficult. e4 is similar to a p-state which is oriented along 
the [001] direction. The electron state e$ is d-like. 

Figure|fjlshows the isosurface plots of the squared wave- 
function |<!>i(r)| 2 for the lowest five hole states h\-h^ with 
and without strain. The light and dark isosurface levels 
are again selected as 0.1 and 0.5 of the maximum prob- 
ability density, respectively. Our atomistic calculation 
shows that the hole states cannot be classified by s-like 
(hi), p-like (/12 and /13) or d-like (/15) shape according 
to their nodal structures. With and without strain the 
hole states underly a strong band mixing. So the cal- 
culated hole states show no nodal structures. Therefore 
the assumption of a single heavy-hole valence-band for 
the description of the bound hole states in a EQD even 
qualitatively yields incorrect results. In contrast to quan- 
tum well systems, the light-hole and heavy-hole bands are 
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Electron states for the b — 10 a base-length dot 

With Strain 



e\ e2 e-i e± es 




Without Strain 

ei e 2 e-i e 5 



FIG. 5: Isosurfaces of the squared electron wavefunctions with and without strain for the embedded b = 10 a pyramidal QD. 
The light and dark surfaces correspond to 0.1 and 0.5 of the maximum probability density, respectively. 

Hole states for the b = 10 a base-length dot 

With Strain 



hi hi /13 hi hs 




Without Strain 



hi hi hz hi h 5 




FIG. 6: Isosurfaces plots of the squared hole wave functions with and without strain for the embedded b — 10 a pyramidal QD. 
See caption of Fig. [S] for more details. 



strongly mixed in a EQD. This result is in good agree- 
ment with other multiband approachea 23 i 24 i 25 i 26 i 27 i 28 i 29 . 

From Fig.Elwe can also estimate the influence of strain 
on the different hole states. Without strain the states hi 
and hi are only slightly elongated along the [110] and 
[110] direction, respectively. Due to strain these states 
are clearly elongated along these directions. The states 
/13-/15 are only slightly affected by strain. 

Another interesting result is that strain effects shift 
the electron states to lower energies and the hole states 
to higher energies as displayed in Fig. 0] Figure 0] also 
reveals that the WL ground-state for electrons and holes 
is shifted in a similar way due to strain. We observe that 
strain decreases the EQD gap E^P = E ei — Ef ll by about 
1.4%, lowering it from the strain-unaffected value 2.12 eV 
to the value 2.09 eV. For a biaxial compressive strain in 
a zinc blende structure, the conduction-band minimum 
of a bulk material is shifted to higher energies while the 



energy shift of the valence-band maximum depends on 
the magnitude of the hydrostatic and shear deformation 
energies^. So one would expect that the electron states 
are shifted to higher energies due to the fact that CdSe 
is compressively strained in the ZnSe-Matrix. This is in 
contradiction to the behavior we observe here. To inves- 
tigate the influence of the WL states on the one-particle 
spectrum we use the same model geometry as shown in 
Fig. |21 but with a considerably smaller WL thickness of 
only one monolayer (ML). A 1 ML thick WL was also 
used before by Santoprete et alw^, Stier et alw& and 
Wang et td~ for an InAs/GaAs EQD. Figure [7| shows 
the comparison of the results for a strain-unaffected and 
a strained pyramidal CdSe EQD with a 1 ML thick WL 
and a base length of b = 10 a. 

On the left-hand side of Fig. [3 the first five electron 
and hole-state energies for an unstrained EQD are dis- 
played while the right-hand side shows the energies for 
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2.8201 =ZnSe CBM 



0=ZnSe VBM 

Nostrain 



WL 


la 


1 ML 




No Strain Strain No Strain Strain 


A eiiE2 [meV] 


162.8 


161.5 


204.1 


221.2 


A e2 ,e 3 [meV] 


0.43 


7.1 


0.5 


10.9 


A hlth2 [meV] 


5.76 


3.7 


7.25 


7.66 


A h 2 ,h 3 [meV] 


16.36 


12.3 


19.66 


15.11 


[eV] 


2.12 


2.09 


2.21 


2.21 



TABLE III: Energy splittings for electron and hole bound- 
states in case of different WL thicknesses. The influence of 
strain-effects on the splittings A ei ,e 2 = \E ei — E e2 |, A e2l e 3 = 
\E e2 — E E3 1, Ah 1 ,h 2 = \Eh t — Eh 2 1 and Aj, 2i j, 3 = \Eh 2 — Eh 3 \ 
is also displayed. The WL thickness is 1 a and one monolayer 
(1ML), respectively. The base length of the pyramid is b = 
10 a. 



Strain 



FIG. 7: First five electron and hole state energies for the em- 
bedded pyramidal CdSe QD with b = 10 a and a one mono- 
layer (1 ML) thick WL. On the left-hand side the results for 
the unstrained EQD are shown while on the right-hand side 
the results for the strained EQD are displayed. The zero 
of the energy scale is the bulk ZnSe valence-band maximum 
(VBM). The energies are compared with the conduction-band 
minimum (CBM) of the bulk ZnSe. 



the strained EQD. For a 1 ML thick WL the lowest elec- 
tron state is, by strain effects, shifted to higher ener- 
gies. This is what one would expect for biaxial com- 
pression of the bulk material. Furthermore the splitting 
of the p-like states e<i and e% is larger compared to the 
results for a 1 a thick WL. The splitting Ag 2 of the 
unstrained EQD with a la thick WL is A° 2 _ e3 = 0.43 
meV whereas for the system with a 1 ML thick WL 
one has A° = 0.5 meV. So the spitting A° is in- 
creased by about 16%. With strain-effects, the splitting 
for the system with 1 ML WL thickness Af 2 ™ n = 10.9 
meV is about 54% larger than the splitting in the sys- 
tem with la WL thickness Af 2 ra ™ = 7.1 meV. Also the 
energy splitting A eii62 between the ground state e\ and 
the first excited state e<i is strongly influenced by the WL 
thickness, namely A ei:E2 = 162.8 meV for the unstrained 
sytem with la WL but A eiiE2 = 204.1 meV for a 1 ML 
WL; with strain effects the splitting A eijG2 is increased 
by about 27% if the WL thickness is decreased from 1 a 
to 1 ML. The results are summerized in Table ITTTl This 
effect mainly arises from the fact, that the bound states 
inside the dot are also coupled to the WL-states. For 
a 1 a WL the wave functions of the bound states show 
also a probability density inside the WL. For a thinner 
WL the leaking of the states into the region of the WL 
is much less pronounced. In this case, the microscopic 
structure inside the EQD and also the strain-affect are 
much more important. This explains the larger energy 
splittings in case of the 1 ML WL compared to the results 



for a la WL. The hole states are influenced in a similar 
manner. In case of a 1 ML WL the energy spectrum of 
the hole states is shifted to higher energies due to the 
strain-effects. This behavior is similar to the behavior 
obtained from the calculations for a la WL (Fig. E}. In 
the 1 ML WL system the energy splittings A/ ll) / l2 and 
Ah 2 .h 3 for the first three hole-states are larger than the 
values we obtain for the system with 1 a WL. These split- 
tings are also summarized in Table ITTTl The WL thick- 
ness also influences the EQD energy gap E^jP. For a 1 
ML WL the electron-states are shifted to higher energies 
in contrast to the behavior of the hole states (compare 
Figs. 01 and EJl- In case of the 1 ML WL the gap energy 
E^p is only slightly affected by the strain. We observe 
here that the strain has opposite effect for electrons and 
hole states: electron states become shallower, approach- 
ing the conduction-band edge, while the hole states be- 
come deeper, moving away from the valence-band edge. 
The knowledge of the single-particle wave functions 
makes the examination of many-particle effects in EQDs 
possible. The single particle wave functions can be used 
for the calculation of Coulomb- and dipole- matrix ele- 
ments as input parameters. For example the investiga- 
tion of multi-exciton emission spectraSi, carrier capture 
and relaxation in semiconductor quantum dot lasers^ or 
a quantum kinetic description of carrier-phonon interac- 
tions' 1 '' is possible. 



IV. RESULTS FOR CdSe NANOCRYSTALS 

A. Geometry and Strain 

In this section we investigate the single particle states 
of CdSe nanocrystals within our TB-model. These nano- 
structures are chemically synthesizedii*i£. The nanocrys- 
tals are nearly spherical in shapei2iiMi§£ and the surface 
is passivated by organic ligands. Due to the flexible 
surrounding matrix, these nanostructures are nearly un- 
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strained^. The size of these nanostructrues is in between 
10 and 40 A in radiua 11 ! 64 ! 66 ! 67 . 

We model such a chemically synthesized NC as an un- 
strained, spherical crystallite with perfect surface passi- 
vation. The zincblende structure is assumed for the CdSe 
nanocrystal. We neglect surface reconstructions2i*22i£i 
and that the surface coverage with ligands is often not 
perfect, though these effects can be important espe- 
cially for very small NCs. However, we concentrate on 
considerably larger NCs than in the before mentioned ref- 
er enceses. Therefore, unlike previous TB work we con- 
centrate here on size and the size dependence of the re- 
sults obtained for the electronic structure of the NCs. 
The TB-parameters, which describe the coupling be- 
tween the dot material and the ligand molecules, are cho- 
sen to be zero. This corresponds to an infinite potential 
barrier at the surface and is commonly used because of 
the larger band gap of the surrounding material 6 —. An 
alternative approach to treat the ligand molecules is dis- 
cussed by Sapra et al. in Ref. The influence of the 
organic ligands on the electronic structure can also be 
investigated more realistically in the framework of mi- 
croscopic descriptions^SiffiiSi. 



B. Single particle states and comparsion with 
experimental results 

We have performed TB-calculations for finite, spher- 
ical, unstrained NCs of diameter between 1.82 nm and 
4.85 nm (corresponding to 3-8 a, when a » 6.07 A is 
the CdSe lattice constant of the conventional unit cell). 
The finite matrix diagonalizations yield both, the discrete 
eigenenergies and the eigenstates. For the largest NCs (of 
diameter 4.85 nm) results for the five lowest lying electron 
and hole eigenstates are shown in Fig. [SJagain in the form 
of an isosurface plot. The lowest lying electronic state e\ 
obviously has spherical symmetry and can be classified 
as a ls-state. Correspondingly the second state e2 has 
the form of a 2s-state and the states e3 t 4 are p-states and 
es is a cZ-like state. Despite the spherical symmetry of the 
system this simple classification is no longer possible for 
the hole states, however. Even the lowest lying hole state 
h\ has no full rotational invariance, i.e. strictly speaking 
it cannot be classified as being an s-state. This is due to 
the intermixing of different atomic TB-valence electron 
states in the NC. Similarly the higher hole states hi — h± 
cannot clearly be classified as an s- or p-like state. This 
is an effect, which simple effective mass models cannot 
account for, but which will have implications in the cal- 
culation of matrix elements between these states, which 
enter selection rules for optical transitions etc. 

In the case of an ideal zinc blende structure as con- 
sidered here we do not obtain any indications of quasi- 
metallic behavior, i.e. of a non-vanishing (quasi continu- 
ous) spectrum of states at the Fermi energy in contrast 
to previous work (assuming an ideal wurtzite structure 
for CdSe nanocrystals)i2i2S^. This is probably due to 



the fact that this quasi-metallic behavior is due to sur- 
face states in the case when no passivation and surface 
reconstruction is taken into account. These surface states 
are formed by the dangling bonds of unsaturated Se at 
the NC surface, which cause s-states in the band gap 
region 38 . In our simplified and restricted TB s c p 3 -basis 
set these s-orbitals at the anions (Se) are not taken into 
account. Therefore, these surface states, which in real- 
ity and in more realistic models are removed (i.e. ener- 
getically drawn down and filled) due to passivation and 
surface reconstruction, do not occur. 

The discrete electronic states of semiconductor NCs 
are experimentally accessible by scanning tunneling mi- 
croscopy (STM)Si4SSiSi. The tunnel current / between 
the metallic tip of the STM and the CdSe nanocrys- 
tal, which is e.g. epitaxially electrodeposited onto a 
template-stripped gold film, is measured as a function 
of the bias voltage V. The conductance (dl/dV) is re- 
lated to the local tunneling density of states. In the 
dl/dV versus V diagram, several discrete peaks can be 
observed. These peaks correspond to the addition en- 
ergies (charging energies) of holes and electrons. The 
spacing between the various peaks can be attributed to 
the Coulomb charging (addition spectrum) and/or charge 
transfer into higher energy levels (excitation spectrum). 
From these measurements the energy gap as well 

as the splitting A ei)62 between electron ground state e\ 
and the first excited state ei can be determined. 

Alperson et alSL investigated CdSe nanocrystals with 
an STM. Here we compare our calculated energy gap 
i£gap°> which is given by the difference between the elec- 
tron, ei, and hole, hi, groundstate, with measured data 
from Ref. lU Figure ED displays the results for CdSe 
NCs with diameters in between 1.82 nm and 4.85 nm. 
Alperson et alSL compare the STM results (dashed dot- 
ted line) with optical spectroscopy measurements (dot- 
ted line) from Ekimov et ali^i. The overall agreement 
with the TB results is very good, especially for the larger 
NCs. Deviations in the case of the small 2 nm NC arise 
from surface reconstructionsi2i2iiSi which are neglected 
here. When the same calculation is done without spin- 
orbit coupling (TB-NO SO), the energy gap E™£° is al- 
ways strongly overestimated by the TB-model, in par- 
ticular for smaller nanocrystals. So the spin-orbit cou- 
pling is important for a satisfactory reproduction of the 
experimental results. For the calculations without spin- 
orbit coupling, the TB-parameters are re-optimized to 
reproduce the characteristic properties (band gap, effec- 
tive masses) of the bulk material. The re-optimized pa- 
rameters are given in Tabled Taking into account the 
electron spin, the lowest electron state e\ is twofold de- 
generated and s-like. This is consistent with the experi- 
mentally observed doubled in the dl/dV characteristic. 
The next excited level is (quasi) sixfold degenerated. The 
spin-orbit coupling splits this into one twofold and one 
fourfold degenerate stated. In the STM measurement 
Alperson et al& observed such a higher multiplicity of 
the second group of peaks. This behavior has also exper- 
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FIG. 8: Isosurfaces (at 30 % of the maximum probability density) of the squared electron and hole wavefunctions of spherical 
CdSe nanocrystals of diameter d = 4.85 nm for the five lowest states 
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FIG. 9: Energy gap Eg^° as a function of the nanocrystal 
diameter d. Compared are the results from our TB-model 
with (TB) and without (TB-NO SO) spin-orbit coupling, a 
STM (STM)££ and an optical measurement (Optical)^. 




2 3 4 5 

Diameter d [nm] 

FIG. 10: Electron (a) and hole (b) energies as a function of 
the nanocrystal diameter d. For electrons (ei-es) and holes 
(hi-hs) the first five eigenvalues are displayed. Each state is 
twofold degenerated. 



imentally 64 and theoretically 34 been observed for InAs 
nanocrystals. The electron energy spectrum for NCs of 
different diameter is shown in Fig.^J|(a). Here the first 
five electron states e\ — e§ are displayed. Note that every 
state is twofold degenerated due to the spin. 

For the hole states the situation is more complicated. 
Alperson et al& observed a high density of states at 
negative bias. The distinction between addition and ex- 
citation peaks is difficult, due to the large number of 
possibilities and the close proximity between the charg- 
ing energy and the level spacing. For the holes we obtain 



that the first two states (hi, /12) and (/13, /14) are fourfold 
degenerated. The energy splitting of these states is also 
very small. These results are consistent to the observa- 
tions of Alperson et al&. Figure (b) shows the hole 
energy versus diamater d for the spherical CdSe NCs. 
Obviously, for all diameters displayed the states hi — /14 
are almost degenerate, i.e. including spin there is almost 
an 8-fold degeneracy of these states. 
Furthermore the calculated splitting A eiie , 2 = E e . 2 —E ei 
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FIG. 11: Splitting A ei ,e 2 = E e2 — E ei between the lowest 
two electronic states as a function of the nanocrystal diame- 
ter d. The results from our TB-model, with (TB) and without 
(TB-NO SO) spin-orbit coupling, and from an STM measure- 
ment (STM, Ref. Wki are displayed. Besides this results from 
infrared spectroscopy (IR, Ref. 72) and optical methods (Op- 
tical, Ref. l67l) are shown. 



between the first two electron states ei and ei is com- 
pared with experimentally observed results for this quan- 
tity. FieurelTTlshows A ei .„ as a function of the nanocrys- 
tal diameter d. The influence of the spin-orbit coupling 
on our results is also investigated. We have done the 
calculations without (TB-NO SO) and with spin orbit- 
coupling (TB). The results of our TB-model for the split- 
ting A eijE2 are compared with results obtained by STM 67 
and by optical methods (optical )Zi. This splitting A eii e 2 
was independently determined experimentally by Guyot- 
Sionnest and Hines^ using infrared spectroscopy (IR). 
Without spin-orbit coupling the TB-model always over- 
estimates the splitting A ei)e2 . Especially for smaller 
nanocrystals the spin-orbit coupling is very important to 
describe the electronic structure. With spin-orbit cou- 
pling the results of the TB-model show good agreement 
with the experimentally observed results. 

V. CONCLUSION 

We have applied an empirical s c p^ TB-model to the 
calculation of the electronic properties of II- VI semicon- 
ductor EQDs and NCs. Assuming a zinc blende lat- 
tice and (per spin direction) one s-like orbital at the 
cation sites and three p-orbitals at the anion sites, the 
TB-parameters for different materials (here ZnSe and 
CdSe) are determined so that the most essential proper- 
ties (band gap, effective masses etc.) of the known band 
structure of the (3 dimensional) bulk materials are well 
reproduced by the TB band structure. Then a CdSe QD 
(on top of a two-dimensional, a few atomic layers thick 
WL) embedded within a ZnSe matrix is modelled by us- 
ing the TB-parameters of the dot material for those sites 



occupied by CdSe and the ZnSe TB matrix elements for 
the remaining sites; suitable averages have to be chosen 
for intersite matrix elements over and for on-site matrix 
elements on anion (Se) sites at interfaces between QD and 
barrier material. Spherical CdSe NCs can be modelled 
similarly by setting the intersite matrix elements between 
surface atoms and atoms in the monolayer of surfactant 
material to zero. The effects of the spin-orbit interaction, 
the band offsets and for the EQDs strain effects are taken 
into account. 

For the EQD systems the numerical diagonalization 
yields a discrete spectrum of bound electron and hole 
states localized in the region of the EQD. Energetically 
these discrete states are below the continuum of the WL 
states. We have investigated the dependence on the EQD 
size and find that the number of the bound states and 
their binding energy increases with increasing dot size, 
therefore the effective band gap decreases. We have also 
investigated the dependence of the bound eigenenergies 
and their degeneracy on strain and on the thickness of 
the WL. Looking at the states themselves one sees that 
conduction band (electron) states can be roughly clas- 
sified as s-like, p-like, etc. states but the valence band 
(hole) states cannot be classified according to such simple 
(s,p,d) symmetries because they arc determined by a mix- 
ing between the different (anion) p-states. This cannot 
be accounted for by simple effective mass models but it 
will be important for instance for the calculation of dipole 
matrix elements between electron and hole states which 
determine the selection rules for optical transitions. For 
the NCs the whole spectrum is discrete, but in spite of 
the spherical symmetry the hole states do not have the 
simple s,p,G?-symmetry but are intermixtures of atomic 
p-orbitals. Even the lowest hole state has no spherical 
s-symmetry but it is 4-fold (8-fold including spin) de- 
generate. The spin-orbit interaction is very important. 
Including the spin-orbit interaction we obtained nearly 
perfect agreement with experimental results obtained by 
STM for the dependence of the band gap and of the split- 
ting of the lowest electronic states on the diameter of the 
NC. 

Compared to (two-band) effective massi£*2£ and multi- 
band k ■ p-modela 23 i 24 i 25 for EQDs our TB model clearly 
has the advantage of a microscopic, atomistic descrip- 
tion. Different atoms and constituents of the nanostruc- 
ture and their actual positions are considered, and this 
may lead to a reduction of symmetries (for instance the 
Ci v symmetry instead of a C^-symmetry) . This may au- 
tomatically lift certain degeneracies and lead automati- 
cally to a splitting, for instance between e% and e% states, 
whereas an 8-Band-/c ■ p model still yields degenerate &2 
and e-i states 2 ^. The effects of inhomogeneous strain can 
be easily incorporated into a TB model by considering the 
deviations of the actual atomic positions from the ideal 
position in the bulk crystal. Only the (empirical) pseu- 
dopotential treatmen-feSSiSL 2 ^ 2 ^ may be still superior and 
more accurate than the TB approach, but in a pseudopo- 
tential descripion a variation of the wave functions within 
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the individual atoms is accounted for and a large number 
of basis states is required. Therefore, a TB description 
is simpler and quicker and allows for the investigation of 
larger nanostructures without loosing information on the 
essential, microscopic details of the structure. Compared 
to other TB models of QD structures, we do not consider 
free standing, isolated QDs (as in Ref. 1421) but we can de- 
scribe realistic QDs (with a WL) embedded into another 
barrier material. We show here that a reduced s c p^-basis 
is already sufficient for a satisfying reproduction of prop- 
erties like optical gaps, energy splittings, etc., and their 
size dependence. Much larger basis sets, namely a sp 3 s*- 
basis^ or even a sp 3 d 5 s* -b&sis*& were used in previous 
TB-models of EQDs. Our reduced, smaller basis set, of 
course, leads to computational simplifications and allows 
for the treatment of larger QDs. Furthermore, we apply 
our TB-model to different materials than investigated 
previously, namely II- VI CdSe nanostructures, and we 
investigate also NCs, for which excellent agreement with 
experimental STM-results could be demonstrated. 

In the future further applications of our TB-model for 
embedded semiconductor QDs and NCs are planned. Of 
course applications to QDs of other materials, for in- 
stance nitride systems, and other (e.g. wurtzite) crys- 



tal structures are possible. Furthermore, EQDs of other 
shape and size (dome-shaped, lens-shaped, truncated 
cones, etc.) or two coupled QDs or freestanding (capped 
and uncapped) QDs can be investigated. A combination 
with ab-initio calculations is also possible by determining 
the TB-parameters from a first-principles band structure 
calculation of the bulk material. Furthermore the influ- 
ence of surface reconstructions and the surfactant mate- 
rial on the results for NCs should be investigated. Espe- 
cially for small NCs these effects are important. Finally 
matrix elements of certain observables like dipole matrix 
elements between the calculated QD electron and hole 
states can be determined, which are important for selec- 
tion rules and the optical properties of these systems. 
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